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We give a complete analysis of indirect determinations of the top 
quark mass in the Standard Model by introducing a systematic 
procedure to identify observables that receive quantum correc¬ 
tions enhanced by powers of We discuss how to use flavour 
physics as a tool to extract the top quark mass. Although present 
data give only a poor determination, we show how future theoret¬ 
ical and experimental progress in flavour physics can lead to an 
accuracy in Mf well below 2 GeV. We revisit determinations of Mt 
from electroweak data, showing how an improved measurement of 
the W mass leads to an accuracy at the level of 1 GeV. 
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1 Introduction 


The top quark mass (M^) is a key input parameter of the Standard Model (SM). Since the 
top quark is the heaviest particle in the SM, its Yukawa coupling yt is sizeable and plays 
a crucial role in determining the predictions of the theory at the quantum level. A precise 
determination of Mt is crucial for: 

• Stability of the electroweak vacuum. Assuming that no new physics modifies the 
short-distance behaviour of the SM, top-quark loops destabilise the Higgs potential 
creating a deeper minimum at large field value. The measured SM parameters lie 
so close to the critical condition for the formation of the large-held minimum that 
the instability scale can huctuate from 10^° GeV to the Planck scale with a variation 
of Mt of merely 2 GeV [1,2]. Any such small change in Mt can have a substantial 
effect in the evolution of the universe at the inhationary epoch [3] and determine the 
viability of scenarios of Higgs inhation [4]. A more precise determination of Mt will 
add important information to our knowledge of particle physics and cosmology. 

• Supersymmetric predictions for the Higgs mass. Within the Minimal Supersymmetric 
Standard Model, the soft-breaking scale that reproduces the observed Higgs mass has 
a strong dependence on Mt. For tan /3 = 1, the supersymmetry-breaking scale is large 
and roughly coincides with the SM stability scale discussed above. For tan (3 = 20, 
maximal stop mixing and degenerate sparticles, precision computations [5] hnd that 
the supersymmetry-breaking scale varies from 1.7 to 2.5 TeV when Mt is varied by 
one standard deviation around its present best-ht value. 

The most precise quoted value of the top-quark pole mass comes from the combination of 
LHG and Tevatron measurements [6] 

(Mi)poie = 173.34 ±0.76 GeV. (1) 

A theoretical concern about the extraction of Mt from data is that the pole top mass is not 
a physical observable. This means that its experimental determination is done through 
the measurement of other physical observables (hnal-state invariant masses, kinematic 
distributions, total rates) that are especially sensitive to Mt- These measurements are 
compared to the results of theoretical calculations, which are expressed in terms of Mt in a 
well-defined renormalisation scheme. In the context of hadron colliders, the extraction of 
Mt suffers from a variety of effects linked to hadronization that are not fully accountable by 
perturbative QGD calculations, like bound-state effects of the ti pairs, parton showering, 
and other non-perturbative corrections (see [7] for a thorough discussion). In practice, the 
extraction of Mt relies on modelling based on Monte-Garlo generators, and this is why [8] 
refers to Mt in eq. (1) as “Monte-Garlo mass”. Its relation to any short-distance definition 
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of the top mass has an inherent ambignity dne to infrared non-pertnrbative effects, which 
probably amonnt to abont 0.3 GeV. Mnch work is ongoing both on the experimental and 
the theoretical sides to control the size of the errors at this level. 

Alternative methods to extract Mt have been proposed, with the aim of hnding ob¬ 
servables whose prediction is theoretically more robust. One interesting possibility is to 
identify observables that can be computed in QCD beyond the leading order in terms of 
the running top mass evaluated at a sufficiently high-energy scale, so that the perturbative 
expansion is completely reliable. The running top mass is then translated into the pole 
mass by means of a relation now known at four-loops in QCD [9]. This programme has 
been applied to the total inclusive tt cross section [10], from which it was possible to extract 
the following values of the pole top mass: 



172.9 ±2.6 GeV ATLAS [11] 
176.7 ± 2.9 GeV CMS [12] 


( 2 ) 


Although the result is theoretically more transparent, the uncertainties in eq. (2) are still 
significantly larger than that in eq. (1). 

These considerations justify the search for alternative strategies to determine Mt, and 
this will be the subject of our paper. Given that the top is the only quark associated to a 
sizeable Yukawa coupling, loop effects in the SM are potentially very sensitive to Mt. Our 
goal is to identify all processes that receive quantum corrections enhanced by powers of 
Mt (in the limit Mt ^ Mw) and infer Mt from their measurements. 

With the experimental confirmation that the Cabibbo-Kobayashi-Maskawa (CKM) ma¬ 
trix gives an overall successful explanation of the transitions among different quark gener¬ 
ations, the main interest of flavour physics has turned towards the search for new effects 
beyond the SM. Indeed, flavour physics provides a unique tool to explore indirectly new 
physics, in a way often complementary to high-energy probes at colliders. However, in 
this paper we want to argue that new developments are guiding us towards a novel use 
of flavour physics data. On the experimental side, the lack of anomalous signals from 
the LHG suggests that new physics may lie at energy scales much higher than previously 
expected. On the theoretical side, present or upcoming improved calculations of flavour 
processes in the SM are opening new frontiers in precision measurements. In light of these 
developments, in this paper we propose to use the comparison between experimental data 
and theoretical predictions of flavour processes as a way to extract the top quark mass, 
under the assumption that the SM is valid up to very short distance scales.^ 

Our strategy is not new: the history of predicting quark masses from loop-induced 
flavour processes is glorious, with some of these predictions made even before the actual 

^For an earlier attempt to determine the top mass from B-B and the rare kaon decays K'^ — > Miyv, 
Kl ^ tt^vD, see ref. [13]. 


3 



discovery of the corresponding particle. This is the case of the charm-qnark mass, whose 
valne was inferred from theoretical considerations on K-K mixing [14] or of the top-quark 
mass, extracted using B-B data [15]. The use of flavour data for an indirect determination 
of Mt is fairly robust from the theoretical point of view, since it relies on controllable SM 
calculations, in which non-perturbative effects are restricted to a few well-known hadronic 
parameters, now under careful scrutiny by lattice calculations. In this paper, we describe 
the status of the extraction of the top mass from the £t of flavour data, finding (Mjflavour = 
(173.4±7.8) GeV. The uncertainty of this extraction is too large to be competitive with the 
direct measurements. However, taking into account foreseeable progress in perturbative 
and lattice calculations, on one side, and experimental measurements, on the other side, 
our projection for the future is that the error can be brought to about 1.7 GeV. 

In our analysis we use the pole top mass Mt as the physical quantity extracted from 
the hts, deriving it, whenever is needed, from the running MS top mass mt{mt) through 
the 0{ag) perturbative expression given in section 3. This choice is dictated mostly by our 
desire to make the results more transparent and to adopt the same variable currently used 
by experimentalists. However, given that the pole mass, unlike the running MS, suffers 
from an (9 (Aqcd) inherent ambiguity, it may become more appropriate in the future, when 
higher accuracy is reached, to modify this choice, abandoning Mt in favour of fnt{mt). 

An important byproduct of our analysis is that the top-mass extraction can be regarded 
as a well-dehned motivation for improved experimental measurements and theoretical cal¬ 
culations in flavour physics. While the exploration for new-physics effects remains the 
most exciting part of the flavour physics programme, the extraction of Mt dehnes a clear 
and concrete benchmark that can be used to determine the goals that experimental and 
theoretical improvements should aim for. 

With the aim of an indirect determination of the top mass, in this paper we also 
reconsider global hts of electroweak observables, hnding {Mt)Ew = (177.0 ± 2.6) GeV, in 
good agreement with previous studies [16,17]. We hud that the determination of Mt is 
dominated by the measurement of M\y. A reduction of the error in the measurement of 
Mw to about 8 MeV, as foreseeable at the LHG [16], can bring down the uncertainty on 
Mt to 1.2 GeV. 

Most of our considerations would be superseded by a futuristic e+e“ collider operating 
at the ti threshold. Such a collider would allow for an unprecedented determination of 
the top mass. Scans of the ti pair production would reach a statistical accuracy on the 
mass measurement of about 20-30 MeV [18]. Recent N^LO calculations can relate such 
measurements to a well-dehned Mt, with a theoretical uncertainty below about 50 MeV [19]. 

Our paper is organised as follows. In section 2 we present a systematic procedure 
to identify observables sensitive, at the quantum level, to powers of the top mass. We 
discuss present and future top mass determinations from havour data in section 3 and 4, 
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respectively, and from electroweak precision data in section 5. Conclusions are given in 
section 6. 


2 Mt dependence of observables in the heavy-top limit 


The large top Yukawa coupling offers the possibility of reconstructing from SM quantum 
effects. In order to identify the physical observables that are most sensitive to the top mass 
at the one-loop level, we develop here a systematic procedure to extract the leading Mt 
dependence predicted by the SM. We work in the heavy-top limit [20] , in which the masses 
of the W and Z bosons are neglected with respect to Mt. This is achieved by considering 
a gauge-less theory with massive quarks, the Higgs boson h, and 3 Goldstone bosons x 
(related by the equivalence theorem [21] to the longitudinal components of the W and Z), 
where the only quark interaction is 

.^ = yt in + h-c. (3) 


Here yt is the top Yukawa coupling, V is the CKM matrix, and we are working in a basis 
in which both quark mass matrices are simultaneously diagonal. The Higgs doublet II is 
given by 



V 


0 

V + h 


(4) 


where v = 246 GeV is the symmetry breaking scale. We can explicitly write eq. (3) as 


^ = 


yt 


-JX (cos|x|/n) (n + h) tt 


x/2 


+ yt 


sin|x|/n 

\x\lv 


h 

1 H— 

V 


+ {x^tRVti diL + h.c.) 


(5) 


where and are the neutral and charged Goldstones, and |yp = -f 2y+y“. The 
next step is to integrate out the top quark using the interactions in eq. (5). The top-less 
effective theory will contain a set of effective operators whose coefficients readily describe 
the leading top-mass dependence in the large Mt limit. 


Ap 

At the level of dimension-4 operators, the hrst diagram in £g. 1 leads to a wave-function 
renormalisation of the Goldstone helds that violates the custodial SU(2) symmetry under 
which X transforms as a triplet. Simple power counting shows that this correction is 
0{yt so we expect a quadratic sensitivity to Mt. Indeed, explicit calculation of the 
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Figure 1: Feynman diagrams illustrating the effective operators generated by integrating 
out the top quark. Also shown is the power counting estimate of their sensitivity to the top 
mass. Dashed lines denote the Higgs boson (h) or the Goldstones (x); solid lines denote 
the quarks. 


diagram in fig. 1 (together with a one-loop diagram obtained from the x'^tt vertex, needed 
to cancel contributions at zero external momentum) reproduces the well-known result for 
the correction to the parameter p = M^/cos6*^M| = 1 -|- Ap 


?>yl 3 GfM^ 
327r2 8V27r2 ' 


( 6 ) 


Z 66, K —TTW and Bg —)■ 


At the level of dimension-5 operators, the second diagram in £g. 1 leads to an effective cou¬ 
pling di'y^dL^d^x^) between a left-handed down current and the derivative of the neutral 
Goldstone which affects the Z couplings. By power counting we estimate the coefficient 
of the dimension-5 operator to be of order \ Vtd\^yl/Mt), which corresponds again to a 
quadratic sensitivity on M*. Explicit calculation of the diagram in £g. 1 gives a correction 
to the Zdidj vertex 





( 7 ) 
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The coupling is defined from 


cos 0w 


i(gtPL + atPR)Ui, 


where in the SM at tree level 


9l = 1 -d + 


1 sin^ e 


w 


5ij , 


i ■ sin^ 6'w . 

9 k = —^—% 


( 8 ) 

(9) 


The vertex correction in eq. (7) gives a quadratic sensitivity to Mt in the Z ^ bb decay 
width 

T(Z^ bb) = [{gf + Agff + ( 9 “)^] , ( 10 ) 

Try/ 

(here given for simplicity in the limit of vanishing bottom mass and neglecting QCD cor¬ 
rections) and in the contribution to the effective Hamiltonians describing K —ttuu and 
Bs 

/\f,sd 

{sLl^dL){y['l^,vi) + h.c., (11) 

(&l7"sl)( 47A) + h.c. (12) 

The effects of eq.s (11)~(12) in the corresponding branching ratios grow as M^. These 
results agree with the leading Mt term of the known full one-loop calculation in the SM. 


AiTiBg and 

The third diagram in fig. 1 leads to a dimension-6 operator involving four di fields. The 
estimate of the coefficient is {Vt*VtjYyt/{IQn'^Mt), exhibiting quadratic sensitivity to the 
top mass. Computing the diagram in fig. 1, we find the AF = 2 interaction 

'^'af =2 = ^ 256*^2 ^2 {diL'^^dji){diLliidjL) + h.c. (13) 

This gives a contribution to CP-conserving and CP-violating observables in meson-antimeson 
mixing with quadratic sensitivity on Mt, in the heavy-top limit. On the other hand, the 
charm-top one loop contribution to ex has no power sensitivity on Mt, in agreement with 
the full SM result. 

Triple gauge boson vertices and WW scattering 

The diagrams in the bottom row of fig. 1 yield a variety of dimension-5 or dimension-6 
operators involving y, h and derivatives, such as h(9^y)^, y(5^y)^, h^(c?^y)^, and X^idf^xY- 
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The usual power counting shows that they have a quadratic sensitivity on the top mass.^ 
These operators contribute to physical observables in triple gauge boson vertices and WW 
scattering. Experimental sensitivity to these effects is too poor to allow for any significant 
determination of Mt. For this reason, we disregard these processes in our analysis, albeit 
their Mf dependence. 

B Xsi 

With the rules of the heavy-top effective theory, it is also easy to identify processes which 
have no power sensitivity on M^. Such processes lead to poor determinations of Mt because, 
in the large Mt limit, one finds at best logarithmic dependences on the top mass. One 
example is S —?• for which the coefficient of the corresponding dimension-6 operator 

mbSLO’^'^bjiF^i, is estimated to be The lack of power sensitivity on Mt 

is confirmed by the full result [22, 23] which, for Mt in the vicinity of its physical value, 
gives 

For this reason, we will not include B —)■ X^y in our analysis. 


Higgs physics 

The heavy-top effective theory also shows that, at present, Higgs physics is not a useful 
player in the game of extracting Mt- The Higgs decays h —)■ yy,yZ can be induced 
by the operators h{d^x^)dyF^i,, whose coefficients are estimated to be 

e^Ht/{IQ'ir'^Mt) (for the first two, which are dimension-5) and ey^/Mt) (for the third, 
which is dimension-6). This corresponds to the well-known result that the amplitudes for 
h —)■ yy, yZ quickly saturate in the large Mt limit. Indeed, from the full SM result we find, 
for Mt around its physical value. 


r(h —>■ yy) oc 


Mt 


173.34 GeV 


0.037 


r(h —)■ Zy) oc 


Mt 


173.34 GeV 


0.014 


(15) 


For the same reason, also h gg offers negligible sensitivity to variations of Mt around 
its physical value. 

Another potential effect comes from the dimension-5 operator h(c}^x)^, generated by 
the first Feynman diagram in the bottom row of fig. 1, whose coefficient is 0{y'^/IQi^'^Mt). 

^There is also an (b(j/^/167r^) correction to hh —i hh scattering and to the triple Higgs coupling. 
The sensitivity of the Higgs self-coupling to yf at the quantum level explains the importance of the 
top-mass measurement for vacuum stability considerations and for the calculation of the Higgs mass in 
supersymmetry. 






An explicit evaluation of the diagram gives the following correction to the Higgs decay 
width into weak gauge bosons 

AT{h^WW*,ZZ*) ^ 

T{h^WW*,ZZ*) 327r2 ’ ^ ^ 

which agrees at the leading order in yt with the known SM result [24]. Even a futuristic 
measurement of the branching ratio at 1% could not determine Mt with an error better 
than 50 GeV. The decays h ^ ZZ, WW, in spite of their quadratic sensitivity on the top 
mass, in practice give no probe of Mt because they are dominated by tree-level effects. 

The process in which the Higgs is radiated off a ti pair offers a direct measurement 
of the top Yukawa coupling. However, the predicted precision in the determination of 
the ratio between the Higgs couplings to top and gluon is in the range 13-17% for the 
LHC with 300 fb“^ and 6-8% at HL-LHC with 3000 fb“^ [25]. This will never become 
competitive with other methods for extracting Mt available in the future. More interesting 
is the case of a hadron collider at 100 TeV, where studies of the ratio tih/tiZ could lead 
to a determination of the top Yukawa with one-percent accuracy. 

We conclude this section by remarking how our analysis based on the heavy-top effective 
theory, after integrating out the top with interactions given in eq. (5), was useful to identify 
the observables most sensitive to Mt- However, for deriving quantitative results on Mt and 
obtain reliable determinations, we have to turn to the full SM expressions of the relevant 
observables. 


3 Extracting Mt from flavour data 


We start by reviewing the basic relations among CKM matrix elements needed for our 
study. Dehning the four parameters X, A, p, r] as 


\Vus\ 

v\yud\^ + \Vus? 


A^m 

X\Vus 


g — if] = 


Vub 

AA3 ’ 


(17) 


the CKM matrix in the Wolfenstein parametrisation [26] becomes 


V 


-X+A‘^X^{l-g-iy) + 0{X'^) 
\AX^{l-Q-ir]) + ^{g+iy) + 0{X'^) 


Ws Wf)\ 

Vcd Vcs Vcb = 

\Vtd Vts Vtb] 

x + o{M) 

l-^-^{\+A^) + 0{X^) 

-AA2+AA4(1 - g-iy) + 0{X^) 


( 18 ) 


AX^{g—iy) \ 
AA2 + C(A^) 
1-^ + 0(A®)/ 
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A=(p,ri) 



Figure 2: The unitarity triangle. 


in agreement with previous results [27]. Unitarity yields the condition 


^udVub ^tdVtb I ^ g 


(19) 


which can be represented as a triangle in the complex plane, see £g. 2 . The vertex A of 
the triangle is given by 


^cd ^cb L ^ 


i - — Q — ir] 


+ C>(A® 


( 20 ) 


while the lengths of the sides CA and BA, denoted by Rh and Rt, respectively, are given 
by 


R^ = 


Rf = 


\VtdVtl\ 


X\Vc 


cb 


i-^ + o{A 


= a /(1 - p )2 + 7^2 = 




td 


The angles f3 and 7 of the triangle are given by the expressions 

sin(/S) sin( 7 ) 


[1 + 0(A-)] 


( 21 ) 

( 22 ) 


Rh = 


Rt = 


sin (7 + /9)’ sin (7 + /3)’ 

which allow us to write the coordinates of the vertex A as 


Q = Rh cos 7 = 1 — cos /3 , r] = Rh sin 'j = Rt sin f3 , 


or, equivalently. 


1 tan 7 11 1 

- = 1 H- - , - =-1-. 

Q tan/? ’ 7 tan /3 tan7 


(23) 


(24) 


(26) 


When searching for new physics, it is customary to determine the four independent 
CKM parameters from tree-level observables, which are presumed to be well described by 
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the SM, and then nse this determination to predict loop processes, which are expected to 
hide new effects beyond the SM. 

In this paper, we take a different perspective: we are assnming the SM to be exactly 
valid and we are interested in extracting Mt from ffavonr processes. We then fix the fonr 
CKM parameters from the most precise measnrements that do not depend on Mt, even if 
they arise at loop leveh: 


IK,.!, |K*|, 7 . / 3 - 


( 26 ) 


The parameters A and A are related to \ Vus\ and \ Vcb\ in the nsnal way, while the expressions 
of g and r] in terms of 7 and are given in eq.s (20) and (25). With this prescription, 
any element of the CKM matrix in eq. (18) can be expressed in terms of the parameters 
in eq. (26). In particnlar, for onr analysis we will need the following combinations 


\VuVi\ 

\yt.v,i\ 

ReAt 

ReAc 

ImAt 


IK.IIW 


smy 


cb 


sin (7 + /3) 

A^ sin (7 - /3) 
2 sin (7 + (3) 


[1 + 0(A‘)] , 


\y 


cb\ 


1 - 


+ C>(A^ 


sin (7 + /3) 


1 + AM - - 


smy 


cos (3 sin (7 + (3) 


+ C>(A^ 


1 - ^ + O(A^) 


-IK 

-ImAe= |Kf,nK 


I siny sin/3 


sin(y + (3) 


1 +A + o(A'' 


( 27 ) 

( 28 ) 

( 29 ) 

( 30 ) 

( 31 ) 


where A^ = (with i = c,t). 


Since the SM predictions for ffavonr observables are often expressed in terms of the 
rnnning MS top qnark mass mt{mt), it is nsefnl to give here the relation between the pole 
top mass Mt and fnt{mt). Acconnting for QCD corrections only^ we find [9] 


Mt 

fritimt) 


1 + 0.4244 as + 0.8345 a] + 2.375 al + (8.49 ± 0.25) a^ = 1.060302(35), (32) 


^An alternative to our choice of CKM input variables could be to take \Vub\ instead of 7 , since present 
relative errors of these two quantities are comparable, see table 1. We prefer the choice in eq. (26) for 
two reasons. First, \Vub\ induces larger uncertainties in the CKM combinations relevant to our analysis. 
Second, 7 is expected to be determined more precisely than \ Vub\ in the future, see table 1 . Thus, we treat 
\Vub\ as a derived quantity, obtained from |Vkh| = |I 4 s||lcb|sind/sin (7 + d)[l + A^/2 + (P(A^)]. 

^This is appropriate for flavour effects where higher-order electroweak corrections have not yet been 
computed. When electroweak corrections have been computed, the result is expressed in terms of Mt, such 
that no conversion in necessary. 
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where = 0.1088. As remarked in the introduction, our choice to express 

results in terms of Mt follows from standard practice, but has the disadvantage of using a 
quantity that it is affected by non-perturbative uncertainties of order Aqcd- 

We can now proceed the to discuss the extraction of the top mass from various flavour 
processes sensitive to Mt. 


A ms. 


The mass differences of the systems in the SM can be written as [43] 

= ^mB^M^13Bjl^y]BSQ{xt)\VtqVt ^\^, q = d,s, (33) 

where rjB accounts for NLO QCD corrections. The LO loop function So{xt) depends on 
Xt = 2yl /where g 2 is the coupling of the SM gauge group SU{2)l and yt is the top- 
Yukawa coupling, and is given by 


So^xt) 


Axt — llxf + Xt 
4(1 - xty 


~ 2 32 ^ 

2{l-xtY^ ' V 173.34 GeVy 


(34) 


The latter equality shows the sensitivity of AMd^s to the top mass in the proximity of its 
physical value. From eq. (33) we obtain the following value for Atub^ 


ArriB, = 


16.9 ± 1.4 
ps 


261 MeV ) 


Mt 


173.34 GeV 


1.52 


0.0401 ) 


Vb 

0.55 


(35) 


Matching this expression with the measurement of Atub^ reported in table 1, we hnd 


=(179.3 ±9.7) GeV. (36) 

Therefore, the current extraction of Mt from Atub^ is affected by an uncertainty of about 
5%. 


Am^^ 


The SM prediction for AruB^ is 


AniB. = 


0.54 ±0.08 
ps 


BbJb^ 


213 MeV 


Mt 


173.34 GeV 


1.52 


and the corresponding determination of the top mass Mt is 


\v.iVS,\ \( m \ 

0.0088 ) V0.55 ) ’ 


(37) 


(M)i,„,^ = (167.0 ±16.8) GeV, (38) 

with an error at the 9% level. Note that the relevant GKM matrix elements and hadronic 
parameters entering Aitib^ are currently less precisely known than those of Aitib^ (see 
table 1) and this explains the smaller error on Mt in eq. (36) than in eq. (38). 
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Observable 

Now (2015) 

Error 2020 

Error 2025 

|GF|xlO-5(GeV-") 

1.16637(1) 

[28| 

- 


- 


Mw{GeV) 

80.385(15) 

[28] 

8 

[16] 

5 

[29] 

Mz{GeV) 

91.1876(21) 

[28| 

~ 



- 

sin^ 6*w 

0.23116(13) 

[28] 

13 

[16] 

1.3 

[16] 


128.952(13) 

[28] 

- 


- 



0.1184(7) 

[28] 

7 

[16] 

7 

[16] 

mc(mc)(GeV) 

1.279(13) 

[30] 

~ 


- 


m_fc(MeV) 

497.614(24) 

[28] 

- 


- 


mB,{MeV) 

5366.8(2) 

[28] 

- 


- 


mB^(MeV) 

5279.2(2) 

[28] 

- 


- 



0.005292(9) 

[28] 

~ 


- 


Ams^(ps"^) 

0.510(3) 

[31] 

- 


- 



17.757(21) 

[31] 

- 


- 


^lf(ps) 

1.607(10) 

[31] 

~ 


- 


iKsI 

0.2249(9) 

[32] 

6 

[33] 

6 

[33] 

\Vcb\x 10^ 

40.9(11) 

[32] 

4 

[33,34] 

3 

[33,34] 

IKfel X 10^ 

3.81(40) 

[32] 

10 

[33,34] 

8 

[33,34] 

sin 2/3 

0.679(20) 

[31] 

16 

[33,34] 

8 

[33,34] 

7 

(73.2+®;3)° 

[31] 

3° 

[33-35] 

1° 

[33-35] 

B{Bs —)■ X10® 

2.8(7) 

[31] 

3 

[33,34] 

1.3 

[33,34] 

^TT^l'l') X 10^^ 

1 7 q+11.5 
' *^-10.5 

[31] 

0.8 

[33,34] 

0.4 

[33,34] 

bIkl 7rVp)xlO^^ 

— 


2 

[33,34] 

0.3 

[33,34] 

|ex| X 10-3 

2.228(11) 

[28] 

- 


- 


/i^(MeV) 

156.3(9) 

[32] 

6 

[33] 

4 

[33] 

Bk 

0.766(10) 

[32] 

7 

[33] 

4 

[33] 

K, 

0.94(2) 

[36] 

? 


? 


Vb 

0.55(1) 

[37] 

0.5 

[38] 

0.2 

[38] 

Vcc 

1.87(76) 

[39] 

? 


? 


Vet 

0.496(47) 

[40] 

? 


? 


Vtt 

0.5765(65) 

[37] 

30 

[38] 

10 

[38] 

6Pc{X)/P,{X) 

0.408(24) [41,42] 

? 


? 


/s.(MeV) 

226(5) 

[32] 

2 

[33] 

1 

[33] 

Bb. 

1.33(6) 

[32] 

2 

[33] 

0.7 

[33] 

JBs/ fBa 

1.204(16) 

[32] 

10 

[33] 

5 

[33] 

Bbs/Bb^ 

1.03(8) 

[32] 

2 

[33] 

0.5 

[33] 


Table 1: Present values and future uncertainties for the most relevant quantities of our 
analysis. In the predictions for future errors we use the symbol when no significant 
improvement is expected, and the symbol “?” when improvement is expected but difficult 
to quantify. 
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The SM prediction for \eK\ can be written as [43] 


\^k\ 


KeGlfj^rriKM^BK . , 

- — - ImAf 

Qy/2TT‘^AmK 


Re\tr]uSo{xt) + ReA^ {rictSo{xc,Xt) 



(39) 


where we have used ImAt = —ImAc, see eq. (31). The multiplicative factor Ks [36] arises 
from long-distance contributions and the parameters rju, rjct, and rjcc accounts for QCD 
corrections. So far, r]u has been calculated at the NLO while rjct and at the NNLO [39, 
40]. The loop function So{xt) is given in eq. (44) and So{xc,Xt) is [43] 


So{Xc,Xt) = Xc 


, xt Sxt 

Xc 4(1 - Xt) 


3xj log Xt 
4(1 - xt)'^ 


2.24 X 10 


-3 


Mt 


173.34 GeV 


0.13 


(40) 


where Xc = m1{mc )/and mc{mc) is the MS charm-quark mass. Inserting the numerical 
values, we hnd 


\£k\ 

KF3 


(1.56 ±0.23) 


f Mt \ 

V 173.34 GeVy 


± (0.50 ±0.19), 


(41) 


which matches the experimental measurement of \sk\ for 


(Mi)|^^l = (185.5 ±22.2) GeV, (42) 

with a 12% error. 


Bs 

The decay Bg —)■ /i+/i“ has been observed by a combined analysis of GMS and LHGb 
data [45]. Although the experimental error is still quite large, see table 1, much progress 
is expected soon. The SM prediction for BR(i?s —)■ /i^/i“) at leading order is [43] 


BR(Rs —)■ ) = 


ml fl ms. 


T 


H 


lOvr^ sin"^ 


' 1 - 


m. 


m 


-\VtgV:,\^Y^{xt) 


(43) 


Bs 


where Q;em(Vfz) = 128.952(13), see table 1, and To(^i) is the loop function 

1.56 


Yoixt) = ^ 


Xt - 4 


± 


3x/ 


_Xt - 1 (xt - 1)' 


log Xt 


0.96 


Mt 


173.34 GeV 


(44) 


The NLO QGD corrections have been included in [46] and found to be very small when 
using the running MS top mass in Yo{xt). The discovery of Bg —)■ has motivated 

^Recently [44] found Kc = 0.963(14) using the most recent lattice QCD inputs. 
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improved SM calculations and NNLO QCD and NLO electroweak corrections have been 
computed [47]. Updating the numerical result of [47] by making use of the input parameters 
of table 1, we hnd 


BR( 5 , ^ /i+/i-) = ( 3.33 ± 0 . 05 ) x 10 “® Rs, 
where Rta and Rs are 

V 0.1184 y V 173.34 GeVy 

R = ( X( X( \VtsVtl/Vcb\ V Tf, 

" V 226 MeVy V 0 . 0409 y V 0.980 ) 1 . 607 ps ’ 


and \VtsV^^/Vch\ is given in eq. (28). Finally, we hnd 

BR(R, ^ /i V) = (3.33 ± 0.24) x 10“® 




173.34 GeV 


3.06 


(45) 


(46) 

(47) 


(48) 


where the uncertainty comes mostly from Vcb and, to a lesser extent, from fs^- Gomparing 
the experimental result for BR(Rs —)■ quoted in table 1 with eq. (48), we end up 

with the following prediction for Mf 

= (163.8 ± 14.7) GeV, (49) 

which suffers from an uncertainty of about 9%. 


K+ 




The branching ratio for —>■ in the SM can be written as [43] 


BR(iF^ —)■ = K+ 


ImAi 


A'(a:.)) 


ReAc , ^ ReA^ , , 

^ ^ ^ Pc{,Xc) + —^X{xt) 


(50) 


where accounts for the hadronic matrix element, which can be extracted from the 
semi-leptonic decays of Kl and Ks mesons [48], and electromagnetic corrections 


k+ = (5.155 ±0.025) x 10 


-11 


A 


0.2249 


(1 ± Aem) , 


(51) 


with Aem = —0.003. X{xt) and Pc{xc) are the loop functions for the top and charm quark 
contributions. The value of Pc{xc) is given by 

Pc{xc) = P^ixc) + 6Pc,u = 0.408 ± 0.024, (52) 
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where P®°(xc) = 0.368 ±0.013, obtained from the results of [41] using the inputs of table 1, 
and SPc^u = 0.04 ± 0.02 [42] arise from short-distance (NNLO QCD and NLO electroweak 
corrections) and long-distance contributions, respectively. On the other hand, the loop 
function X{xt) can be written as 


(y. (y. 

X{xt) = Xo{xt) ± -^Xx{xt) ± —X^^{xt ), 

47r 47r 


(53) 


where X^^xt) accounts for the LO result [43] 


M^t) = ? 


Xt + 2 3xt-Q 


1.50 




173.34 GeV 


1.15 


(54) 


while Xi{xt) and Xew{xt) are relative to NLO QCD and electroweak corrections, respec¬ 
tively. 

The full two-loop electroweak corrections to the top-quark contribution Xt has been 
computed [49] , bringing the theoretical uncertainty related to electroweak effects well below 
1%. A very accurate approximation of the full result is captured by the expression [49] 


X{xt) = 




Vx 


47r 


Mt 


A — B (7 173.34 GeV -|- D 


M, 


173.34 GeV 


Xo(xt), 


(56) 


where r]x = 0.985 stems from NLO QCD corrections, while A ~ P ~ 1.12, C ~ 1.15, and 
Il~0.18 arise from NLO electroweak corrections. 

Using the inputs of table 1, we hnd the following prediction for BR(iL+ —>■ 7i~^ui>) 


BR(A:+ ^ TT+z/h) = (8.42 ± 0.61) x 10"^^. 


(56) 


Even if —)• vr+z/P has been already observed, its experimental resolution is so poor (see 

table 1) that any extraction of Mt from —)• is meaningless at present. For this 

reason, we postpone the determination of Mt to the next section, where we discuss future 
theoretical and experimental improvements. 


Kl -a- 


The branching ratio for Kl —)• vr'^z/P in the SM is fully dominated by the diagrams with 
internal top exchanges, with the charm contribution well below 1%. It can be written as 
follows [43] 


BR(iLi —)■ 7r°z/z/) = kl 


ImA 


-X{xt) 



(57) 


where kl accounts for the hadronic matrix element and is given by [48] 


Kl = (2.223 ±0.013) x 10"^° 


A 


0.2249 


(58) 
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Due to the absence of the charm contribution in eq. (57), the theoretical uncertainties in 
—)■ 7 r°z/z/) arise only from the CKM matrix elements. We hnd that the current 
prediction for BR( 772 . is 

BR{Kl vrVz/) = (2.64 ± 0.41) x 10“^^, (59) 

where we have used the inputs of table 1 . 

This process has not been observed yet. Future prospects for the extraction of Mt from 
Kl ^ will be addressed in the next section. 


Global fit 


The determinations of Mt from the various flavour processes and their combination are 
summarised in £g. 3. Our result for the pole top mass extracted from flavour physics is 


(M0flavour = (173.4 ±7.8) GeV 


( 60 ) 


This result is compatible with the collider determination in eq. (1), but the error is too 
large to be competitive. 

In principle, the extraction of (Mt)fiavour would require a global £t of all flavour observ¬ 
ables in which the CKM parameters and the top mass are allowed to float independently. 
However, in practice, our procedure of hxing the CKM parameters in eq. (26) from pro¬ 
cesses that are insensitive to Mt and then determine Mt from the remaining observables 
is perfectly adequate and leads to results identical to those from a global £t. Actually, as 
shown in £g. 3, the determination of Mt is dominated by Ams^, which depends on the 
CKM parameters only through the combination Equation (28) shows that this 

combination is equal to \Vcb\, up to a dependence on the angles 7 and (3 suppressed by two 
powers of A. This means that essentially \Vcb\ alone drives the error on the determination 
of Mt attributable to CKM elements, while the less precisely known parameters 7 and (3 
play only a minor role. As we will show in the next section, Bg —)• fMiJi~ will soon become 
an equally important process for the determination of Mt and its CKM dependence, as in 
the case of is given by So our conclusion that \Vcb\ is the most important 

CKM parameter for Mt extraction is likely to hold true even after future theoretical and 
experimental improvements. Let us turn now to discuss our forecast for the future of Mt 
determinations from flavour processes. 


4 Future determinations of Mt from flavour 

The current determination of Mt from flavour processes in eq. (60) will soon improve 
thanks to upcoming experimental and theoretical progress. Figure 4 (left panel) shows 
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Present fits and projections for 2020 (red) and 2025 (green) 


166.6+ 15.9 

AniB^ 

178.7 + 9.8 


163.8 + 14.1 

K* -->^71* vv 

179.2 + 166.6 

185.6 + 22.9 

Global flavor fit 

173.4 + 7.8 


direct 173.3 + 0.8 

_ ^^^^^ ^^^^ ^^^^^^^^^^^^^ _ 

150 160 170 180 190 200 

Pole top mass M, in GeV 

Figure 3: Summary of present and future determinations of Mt from flavour data. For 
future projections, we have fixed the central value of Mt to the present direct measurement. 


how the uncertainty on the value of Mt extracted from the global £t changes, as we vary 
the uncertainties of each observable one at a time. We only show the effect of the input 
parameters that have a signihcant impact. We see that more precise measurements of 
BR(i?s —)■ iMii~) and Wb, and a more precise computation of fs^ are the key elements 
for improvements in the determination of Mt. 

However, future improvements will come simultaneously from many observables. Thus, 
in this section we estimate the future situation, in light of new measurements from LHCb, 
Belle II, and NA62, progress in unquenched lattice QCD calculations, as well as improve¬ 
ments in theoretical calculations of QCD and electroweak short-distance effects. We will 
outline the error budget of each flavour observable aiming to quantify the improvements 
needed to bring the error on Mt at the 1% level. 

Future determinations of Mt from Atub^ 

In the case of AmB„, the error budget is 

5(AmsJ = (±1.07^02^^^ ± 0.91|y^,| ± 0.31,,) ps"', (61) 

where, hereafter, we always assume the current SM central values. Eq. (61) clearly shows 
that the major sources of errors arise from (in order of importance): i) the hadronic 
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M, from flavour data 


M, from electroweak data 



Figure 4: Uncertainty on the value of Mt extracted from flavour (left) and electroweak 
data (right) as the uncertainties on some key observables are varied one at a time. 


parameters ii) the CKM matrix elements |l^s| ~ I Kb I , and in) short-distance 

QCD effects encoded in the parameter rfB- On the other hand, the experimental error on 
ArriB^ is ±0.02exp, thus, totally negligible. 

The goal of lattice QCD, concerning the calculation of quantities related to flavour 
physics, is to reach a resolution at the 1% level (or even slightly better) by 2025, see table 
1. However, at this level of precision, one should also consider small effects such as isospin 
breaking and electromagnetic effects, which are 0[{md—mu)/A qqd] and 0{a), respectively, 
and thus at the 1% level. First lattice studies of isospin breaking and electromagnetic effects 
have been performed in the last years leading to very promising results [50-52] . Moreover, 
lattice calculations of form factors of exclusive semileptonic H-decays are crucial to extract 
I Kb I and I Kb I- They are extracted from more noisy three-point correlators and imply 
an extrapolation in the transfer momentum, which is computationally intense. For the 
semileptonic decays B — D/D*fv, however, one measures on the lattice the difference of 
the form factor from unity (he. the SU(3) or heavy-quark symmetric limit), so that the 
uncertainty on the form factor itself turns out to be smaller. 

Concerning |Kb|, there are discrepancies between its inclusive and exclusive extrapola¬ 
tions from tree-level decays and Belle II should resolve this problem.® Overall, exclusive 
determinations are expected to be more precise because they are easier to perform experi- 

®In our analysis, we use the average of inclusive and exclusive determinations of ref. [32], see table 1, 
which does not include the most recent calculation of | Feb |inci of ref. [53]. 
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mentally and also because the calculations of the relevant form factors from lattice QCD 
are less challenging than in the inclusive case. Finally, the error associated with t]b will be 
reduced significantly, at least by a factor of 3-4 once NNLO QCD and NLO electroweak 
calculations will be available. 

Assuming the expected improvements by about 2025, see table 1, we have 

5(AmsJ = ± 0.25|v^,| ± 0.06,,) ps"', (62) 

which corresponds to about a factor of 4 improvement in the overall error compared to the 
current error, see eq. (61). We can determine imposing that the SM prediction 

for AniB^ matches its experimental measurement. This leads to the relation 


S{Mt) 




Mt 


±0.66 


\ 


Btfo, 


±4 




where SAmsg refers to the experimental uncertainty on Am^,. Finally we hnd 






±3.6 GeV 
±2.1 GeV 


( 2020 ) 

(2025) ’ 


(63) 


(64) 


in good agreement with our numerical results in £g. 4. These values have to be compared 
with the current uncertainty = ±10 GeV, see eq. (36). 


Future determinations of Mt from 

In the case of Am^^, the current error budget is 


5(AmBj= ±0.056^1/2, ± 0.029|v;j ± 0.001^3 ± 0.048, ± 0.010,, ps-Q (65) 

\ -“Sd JBd J 

while the experimental error ±0.003exp is negligible. Many considerations done for Am^^ 
hold here too, the only difference being that the uncertainties on ^ind are 

larger than in the Aitib^ case. Assuming the expected improvements by 2025, see table 1, 
we have 


^(AmBj = \^± 0.008^v 2^^ ± 0.008|v^,| ± O.OOl;? ± 0.007, ± 0.002,, J ps"^, (66) 

which corresponds, as in the ArriB^ case, to about a factor of 4 improvement compared 
to the current uncertainty. Notice that now the experimental error ±0.003exp is no longer 
negligible. The projected errors on (5(±h)^^ by 2020 and 2025 can be found from eq. (63) 
and read 


SiM) ~l ±6.6 GeV (2020) 

V ~ j ^3 J Qj.y po25) ’ 


(67) 
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in good agreement with our numerical results in fig. 4. These values have to be compared 
with the current uncertainty = ±16 GeV, see eq. (38). Therefore, by around 

2025, the expected uncertainty on the value of Mt extracted from will be about 

1 . 6 %. 


Future determinations of Mt from Bg 

In the case of BR(i?s —)■ not only theoretical but especially experimental uncer¬ 

tainties have to be reduced significantly in order to extract the top mass with an improved 
accuracy. On the experimental side, the LHCb collaboration aims at reaching a 10% res¬ 
olution on BR(i ?5 —in a few years. The final goal, after the LHCb upgrade, is a 
resolution around (4 — 5)%. On the theoretical side, the main sources of uncertainties arise 
from the decay constant and \Vcb\- The error budget for BR(i?s —)■ /i^/i“) is 

- ^ ' = ±0.05.k ± 0,15/,. ± 0.18|v,.| ± 0.02,j , (68) 

where ± 0.06th stems from the estimated error from higher-order effects, as discussed in [47]. 
On the other hand, the experimental error ±0.84exp is by far dominant at present. The 
situation is expected to improve greatly in the future. By 2025 the error budget will be 

^ ^ ± O-Olth ± 0.03/^^ ± 0.05|v;,| ± 0.02^1^, (69) 

assuming that the errors from higher-order effects will be significantly reduced. 

Matching the SM prediction, see eq. (45), with the experimental result leads to the 
determination of the top mass uncertainty through the relation 








= ±0.33 


IBs 


±4 




cb\ 


\Vc 


cb\ 


± 


OTh 


T\ 


H 


± 


6B 


fi/i 


B 


ii/i 


(70) 


where 5B^^ stands for the experimental error on = BR(Rs —>■ ). We predict. 


±5.3 GeV (2020) 
±2.4 GeV (2025) ’ 


in good agreement with our numerical results in fig. 4 


(71) 


Future determinations of Mt from ek 

In order to reduce significantly the determination from ex, one would need to improve 
especially the uncertainties on t]cc, f]ct and Vcb- It is important to stress that a final answer 
about the errors on t]cc and r]ct are expected to come from lattice QCD calculations. In 2-3 
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years, a fully controlled calculation reducing the total error coming from rjcc and rjct to the 
1 % level should be available, although this is a challenging task for lattice simulations. 
Let us now study the current error budget of ex which is given by 

^ = ± 0.17|y^,| ± 0.14^ ± 0.05;3 ± ± 0.15,,, ± 0.08„, ± 0.02,,, ± 0.04,,,. (72) 

On the other hand, the expected error budget by 2025 is 

^ = ± 0.05|v^„| ± 0.02, ± 0.02^ ± 0.02^v.^^ ± 0.02„„„„,„ , (73) 


where we have assumed that the non-perturbative uncertainties encoded in rjcc and rjct 
will almost disappear thanks to lattice calculations [54]. We hnd the following top mass 
uncertainties 

r ± 8 GeV (2020) 

I ± 5 GeV (2025) ’ 

in good agreement with our numerical results in fig. 4. 




Future determinations of Mt from 


In the case of BR(iL’'' —>■ 7r~^i>9), the by far dominant uncertainty comes from \Vch\ and 
to a lesser extent from the long-distance effects encoded in Pc{X). Goncerning the latter 
uncertainty, there is ongoing activity by lattice QGD collaborations aiming to reduce it to 
the 1% level in a few years from now. 

On the experimental side, the NA62 experiment at GERN aims to measure BR(iL’'' —)■ 
7r~^uu) with a 10% accuracy by 2018 while a 5% resolution could be the hnal goal of NA62. 
The current error budget for BR(iL’'' —)• vr+i/p) is 


5BR(A:+ ^ Ti+uP) 

RRii 


± 0.52|y„,| ± 0.43, ± 0.02^3 ± 0.23p,. 


By 2025 the error budget will presumably be 


(75) 


5br(a:+ ^ TT+i/p) 


±0.14|y,,l ±0.07,±0.01^, 


(76) 


where we have assumed that the non-perturbative uncertainties encoded in Pc will dis¬ 
appear thanks to lattice calculations [55]. We estimate the following future top mass 
uncertainties 


^(^t)p'+-S.7r+w 


± 12 GeV (2020) 
±7 GeV (2025) ’ 


( 77 ) 


in good agreement with our numerical results in fig. 4. 
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Future determinations of Mt from Kl 

On the experimental side, the KOTO experiment at J-PARC plans to reach the SM level 
for BIi{KL —)■ 7 r°i/p) in a few years from now. The expected data corresponds to a few 
events. With an upgrade of the KOTO experiment the hnal goal is to obtain a sample of 
about 100 SM events, corresponding to a 10% resolution on BIi{KL —)■ 

On the theoretical side, since BIi{Ki —)■ ir^uu) is fully dominated by short-distance 
effects, the main sources of uncertainties arise from iWbl, and 7 , see eq.s (31),(57). In 
particular, the current error budget for BR(iPi —)■ tt'^z/P) is 


6BR{Kl ttVp) 

By 2025 the error budget will be 

6BR{Kl ttVp) 


10 


-11 


= ± 0.28|v;,I ±0.23.^ ±0.19^3 • 


= ± O.OSiy^j^i ± 0.04..^ ± 0.08^ . 


(78) 


(79) 


We estimate that the top mass uncertainty will be 


S{Mt) 




M, 


= ±0.43 WI 6 




ch 


\Vc 


cb 


±4 


5Rt sin (3 
Rt sin (3 


± 


5K, 


(80) 


where is the experimental error on = BR{Kl —)■ Our projection is 


5{Mt) 




±57.5 GeV (2020) 
±9.2 GeV (2025) 


(81) 


Global fit 

The expected determinations of Mt by 2020 and 2025 from the various flavour processes 
and their combination are summarised in £g. 3. From our global £t we predict 

^(Mi)flavour ~ 3.1 GeV (2020), (82) 

5(Mi)flavour ~ 1.7 GeV (2025). (83) 

From £g. 3 we also learn that and Bg —)■ are the most accurate Mt discrimi¬ 
nators, while other observables like ArriB^, and K —)■ play a sub-leading role. The 

latter point is also illustrated by £g. 5, which shows how experimental improvements in 
each flavour observable affect the uncertainty on Mt, assuming present (left) and future 
(right) theory uncertainties. 

We are ready now to summarise the main results of this section. 
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With present theory uncertainties 


With future (s;2025) theory uncertainties 




Fractional experimental uncertainty 


Fractional experimental uncertainty 


Figure 5: Uncertainty on the value of Mt extracted from single flavour observables, as the 
corresponding experimental errors are varied. The arrows mark the present experimental 
uncertainties of the various observables. We assume present theory uncertainties (left) 
and those predicted for around 2025 (right). 


• Since AmB„ and Bg —)■ th® dominant observables and only depend on the 

CKM parameters through the combination |= \Vcb\+0{X^), the determination 
of Mt essentially does not require a complete global £t analysis. 

• A precise determination of Mt from requires substantial improvements of \ Vcb\ 

and see eq. (33). Concerning \Vcb\i a joint effort of experiments and theory 

is necessary. The measurements of i? —)■ D/D*ih> branching ratios by Belle II and 
lattice QCD calculations of the relevant form factors should enable us to extract \Vcb\ 
at or even below the 1% level by around 2025. At the same time B^J^fs^ should 
be calculated by lattice QCD with a precision of about 0.5%, see table 1. At this 
level of precision, it will be mandatory to improve also theoretical calculations by the 
inclusion of NNLO QCD and NLO electroweak short-distance effects. On the other 
hand, the experimental resolution on Am^^, which is already at 0.1%, needs not to 
be improved. As a result, we expect ^{Mt)^^^ ~ ±2.1 GeV by about 2025. 

• Unlike the AmB„ case, a precise determination of Mt from Bg —)■ pApi~ requires 
primarily experimental progress in the measurement of its branching ratio. On the 
theory side, the leading uncertainties stem from fB„, \Vcb\ and, to a lesser extent, 
from higher-order effects which are estimated to induce an error of 1.5% [47]. By 
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= 178.1 ± 3.1 
{M,\^ = 199. ± 15. 

= 286. ± 109. 
(Mt\^ = 158.9 ± 14.0 


(M,)ew = 177.0 ± 2.6 

I—•—I 

(M,)poie = 173.2 ± 0.9 

[-•H 


150 160 170 180 190 200 210 


Pole top mass M, in GeV 

Figure 6 : Summary of present determinations of Mt from electroweak data. 


about 2025, the expected error in BR{Bs —)• driven by the combination of fs^ 

and \Vch\ will be below 2% while the experimental error around 4-5% and therefore 
still dominant. So, we expect ^ ±2.5 GeV by about 2025. In case 

the experimental error on BR(i?s —t should reach the 2% level, it would be 

mandatory to improve the estimated 1.5% error associated with higher-order effects. 
In the latter case, we would obtain ~ ± 1.5 GeV, see £g. 5. 

Even though we have identihed Ams^ and —)■ as the dominant Mt discrim¬ 

inators, we must stress that improvements in the determinations of all other observables 
discussed in this paper are also important. Indeed, our basic assumption for the extraction 
of the top mass from flavour physics relies on the validity of the SM up to large energy. 
In order to establish whether this situation is realised in Nature or not, we need a global 
analysis conhrming that the GKM picture of flavour and GP violation is indeed correct 
also after the expected theoretical and experimental rehnements. 

5 Extracting Mt from electroweak precision data 

Electroweak observables depend on the top mass (and on the Higgs mass) only through 
the £ 1,62 ±3 parameters that describe corrections to the tree-level propagators of the weak 
gauge bosons, and through the Sb parameter that describes corrections to the Zbb ver¬ 
tex [56]. These parameters are related to combinations of physical observables and can be 
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extracted from a global fit of experimental measurements: 


r M 

= +(5.6+ 1.0) X 10-® 


( 1 

0.80 

0.86 

-0.32\ 

1 ^2 
^3 

= -(7.8+ 0.9) X 10-® 

= +(5.6+ 0.9) X 10-® 

with p = 

0.80 

0.86 

1 

0.57 

0.57 

1 

-0.31 

-0.21 

1 Sb 

= -(5.8 + 1.3) X 10-® 


\-0.32 

-0.31 

-0.21 

1 / 


. ( 84 ) 


where p is the correlation matrix/ The SM predictions for these observables, for the central 
values of a^i^Mz) and aem{Mz) and around the measured values of Mt and Mh, are® 

{ £1 = +5.22 X 10-® (Mi/173.34 GeV)®'^® (Mi,/125.09 GeV)-°'^® 

£2 = -7.32 X 10-® (Mi/173.34 GeV)-o-®9 (M;,/125.09 GeV)-°-°® 

£3 = +5.28 X 10-® (Mi/173.34 GeV)-°-°^ (Mi,/125.09 GeV)°-®i ' ^ ’ 

Eb = -6.95 X 10-® (Mi/173.34 GeV)-2-i® 


As discussed in section 2, in the large Mi limit the one-loop corrections to Ei = Ap 
and Eb = —2Ag^l grow as M/, while E 2 and £3 only have a milder In Mi dependence. 
Furthermore £1 and £3 have (in the large Mh limit) a In Mh dependence, which leads to a 
negligible uncertainty, now that Mh = (125.09 + 0.24) GeV is precisely measured. Figure 6 
summarises the various single determinations of Mt from the e pseudo-observables, from 
which we derive our result of the global electroweak fit: 


(Mi)Ew = (177.0 ± 2.6) GeV 


( 86 ) 


This result agrees with recent global fits that found Mt = (177.0 ± 2.4) GeV [16] and 
Mt = (176.6 + 2.5) GeV [17]. In £g. 4 (right panel) we show how the uncertainty on the Mt 
determination from the global £t changes, when uncertainties on the various observables 
are changed one-by-one. We only show the effect of those that have the most significant 
impact. We see that: 


• The measurement of Mw plays the key role, since we hnd 6Mt/Mt = 696 M]y/Mw. 
This means that measuring M^ with a precision of 8 MeV (as foreseeable after 
combination of the full LHG dataset [16]) can lead to a determination of Mt within 
about 1.2 GeV. On the other hand, measurements of the WW production cross 
section at the ILG could reduce the error on Mw to about 5 MeV [29], corresponding 
to a determination of Mt at the level of 0.7 GeV. 

^The mean values gt, the errors at and the correlation matrix pij determine the as 

“ M»)(fr^)/1(+ - +), where (CT^)y = a.pi^aj . 
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• The determination of Mt can also be improved by better measurements of the various 
asymmetries (blue lines in the right panel of £g. 4), which determine the weak mixing 
angle, and of Ri,. However, only with a reduction of the present errors on these quan¬ 
tities by more than a factor of 3 one can start observing meaningful improvements 
on the determination of Mt. 

• The £t is not crucially sensitive to other parameters. In particular, the uncertainty 
on Mt would be affected only if the error on aem{Mz) were underestimated by more 
than a factor of 2 . 


Since the determination of Mt from the global fit of electroweak data is largely domi¬ 
nated by Mw, it is useful to reconsider the extraction of Mt using Mw as the only input 
quantity. The value of Mw enters the definition of the pseudo-observable Ar^y, which is 
defined as the ratio of two different determinations of the weak angle: 


Arvy = 1 — 


/ a/2Gf^I 


= (-25.4 ± 0.95m., ±0.10„„) x 10 


-3 


(87) 


The numerical value has been obtained by taking the experimental values of the SM pa¬ 
rameters given in table 1. Equation (87) shows that the uncertainty in aemiMz) has a 
subdominant influence with respect to Mw- The quantity Arw can be computed in the 
SM and expressed in terms of the £ parameters as 


Arw = 


tan ^ 6 *w + (tan 6 *w — 1 ) + 2 £3 

Mt \ ( Mh '' 


= -24.0 X 10 


-3 


173.34 GeV 


125.09 GeV 


( 88 ) 


The above numerical expression has not been obtained by simply replacing eq. (85) into 
eq. ( 88 ), but rather by using the full two-loop result that can be extracted from the calcu¬ 
lation presented in [57]. Such result has never been included before in global electroweak 
hts. By comparing eq. (87) with eq. ( 88 ) we hnd 


{Mt)M^ = (177.7 ±2.8) GeV, (89) 

which essentially reproduces the result in eq. ( 86 ), derived from the global £t. This shows 
that the determination of the top mass from electroweak data is almost completely driven 
by Arw and a full global £t is superfluous if one is interested in obtaining a simple, but 
reliable, estimate of Mt. 


6 Conclusions 

In this paper we have analysed indirect determinations of the top quark mass Mt. For 
this purpose, in section 2 we have presented a systematic procedure to identify observables 
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that, at the quantum level, have power sensitivity on the top mass, in the limit Mt 3> Mw- 
This is done by considering an effective theory obtained after integrating out the top quark 
in the gauge-less limit of the SM. We have divided the physical quantities sensitive to Mt 
into two classes: flavour observables and electroweak observables. 

In section 3 we have discussed how the top mass Mt is determined through Mt- 
dependent quantum effects in the physical quantities 

ArriB, , Bs -)■ -t vr+z/h, Kl -)■ . (90) 

The determination of Mt from the first two observables essentially requires only Vet as 
CKM input. Moreover, these two observables provide the best probe of Mt among flavour 
processes. Hence, ArriB^ and Bg —)■ combined with a determination of Wfe and the 

lattice parameters amj /s^, are sufficient to extract a fairly accurate estimate of 

the Mt determination from flavour physics. Adding to the analysis the other observables 
listed in (90) requires a complete joint fit with all CKM parameters and has a limited 
impact on the extracted value of Mt. Our results are summarised in fig. 3: at present 
flavour data determine Mt = (173.4 ± 7.8) GeV, with Attib^ and Bg —)■ being the 

best toppometers. 

In section 4 we have discussed how the uncertainty on Mt from flavour determinations 
is expected to decrease significantly in the future, mostly thanks to better measurements 
of Bg —)■ to better lattice computations of the hadronic parameters entering the SM 

prediction of Attib^ and Bg —)■ iJi~ and to improved theoretical calculations of short- 
distance effects. We have estimated that the uncertainty on Mt can be brought down to 
3 GeV by 2020 and to 1.7 GeV by 2025. 

In section 5 we have considered electroweak data, finding that at present they determine 
Mt = (177.0 ± 2.6) GeV. We have found that Mw and V{Z —)■ W) are the most sensitive 
quantities, because of the power dependence on Mt of their quantum corrections. However, 
Mw is by far the best toppometer in electroweak physics. We have presented analytic 
expressions to extract Mt from measurements of Mw which take into account recently 
computed two-loop electroweak quantum corrections [57], not yet included in global fit 
codes. Figure 4 (right panel) shows that a more precise measurement of Mw is the key 
player for an improved determination of Mt from electroweak observables. As experiments 
at the LHG are expected to reduce the uncertainty on Mw to about 8 MeV [16], it is 
foreseeable that electroweak physics will determine Mt with a precision of about 1.2 GeV. 

In the future, a global fit of all indirect determinations of Mt, from both electroweak 
and flavour data, will provide significant information. Even if indirect measurements do 
not surpass direct determinations in precision, the comparison between indirect and direct 
analyses will carry essential information, especially in view of the theoretical ambiguities 
in the extraction of Mt from collider experiments. 


28 


Acknowledgments 

We thank Pier Paolo Giardino for having provided us with the full two-loop prediction for Arw and Satoshi 
Mishima for having provided us with recent computations of the e parameters. We thank P. Gambino 
and T. Mannel for very useful discussions about Vcb, V. Lubicz, G. Martinelli and G. Sachrajda for very 
helpful discussions about lattice QCD, and J. Brod, M. Gorbahn, E. Lunghi and M. Misiak for very useful 
information about short-distance effects. PP thanks G. Buchalla, G. Isidori, U. Nierste and J. Zupan for 
the invitation to the MIAPP workshop Flavour 2015: New Physics at High Energy and High Precision, 
where part of his work was performed. The research of PP is supported by the ERG Advanced Grant No. 
267985 (DaMeSyFla), by the research grant TAsP (Theoretical Astroparticle Physics), and by the Istituto 
Nazionale di Fisica Nucleare (INFN). 


References 

[1] G. Degrassi, S. Di Vita, J. Elias-Miro, J. R. Espinosa, G. F. Giudice, G. Isidori and A. Strumia, 
JHEP 1208 (2012) 098 [arXiv:1205.6497]; D. Buttazzo, G. Degrassi, P. P. Giardino, G. F. Giudice, 
F. Sala, A. Salvio and A. Strumia, JHEP 1312 (2013) 089 [arXiv: 1307.3536]. 

[2] F. Bezrukov, M. Y. Kalmykov, B. A. Kniehl and M. Shaposhnikov, JHEP 1210 (2012) 
140 [arXiv:1205.2893]; A. V. Bednyakov, B. A. Kniehl, A. F. Pikelner and O. L. Veretin, 
[arXiv: 1507.08833]. 

[3] J. R. Espinosa, G. F. Giudice and A. Riotto, JCAP 0805 (2008) 002 [arXiv:0710.2484]; A. Kobakhidze 
and A. Spencer-Smith, Phys. Lett. B 722 (2013) 130 [arXiv:1301.2846]; K. Enqvist, T. Meriniemi 
and S. Nurmi, JCAP 1407 (2014) 025 [arXiv: 1404.3699]; M. Herranen, T. Markkanen, S. Nurmi 
and A. Rajantie, Phys. Rev. Lett. 113 (2014) 21, 211102 [arXiv: 1407.3141]; A. Hook, J. Kearney, 
B. Shakya and K. M. Zurek, JHEP 1501 (2015) 061 [arXiv:1404.5953]; J. Kearney, H. Yoo and 
K. M. Zurek, Phys. Rev. D 91 (2015) 12, 123537 [arXiv: 1503.05193]; J. R. Espinosa, G. F. Giudice, 
E. Morgante, A. Riotto, L. Senatore, A. Strumia and N. Tetradis, [arXiv: 1505.04825]. 

[4] F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B 659 (2008) 703 [arXiv:0710.3755]. 

[5] E. Bagnaschi, G. F. Giudice, P. Slavich and A. Strumia, JHEP 1409 (2014) 092 [arXiv: 1407.4081]; 
J. P. Vega and G. Villadoro, JHEP 1507 (2015) 159 [arXiv:1504.05200]. 

[6] ATLAS and CDF and CMS and DO collaborations, [arXiv:1403.4427]. 

[7] A. Juste, S. Mantry, A. Mitov, A. Penin, P. Skands, E. Varnes, M. Vos and S. Wimpenny, Eur. Phys. J. 
C 74 (2014) 10, 3119 [arXiv: 1310.0799]; K. Kroninger, A. B. Meyer and P. Uwer, [arXiv: 1506.02800]. 

[8] Top Quark Working Group Collaboration, [arXiv: 1311.2028]. 

[9] P. Marquard, A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. Rev. Lett. 114 (2015) 14, 
142002 [arXiv:1502.01030]. 

[10] U. Langenfeld, S. Moch and P. Uwer, Phys. Rev. D 80 (2009) 054009 [arXiv:0906.5273]; V. Ahrens, 
A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, Phys. Lett. B 703 (2011) 135 [arXiv: 1105.5824]; 
M. Beneke, P. Falgari, S. Klein, J. Piclum, C. Schwinn, M. Ubiali and F. Yan, JHEP 1207 
(2012) 194 [arXiv:1206.2454]; M. Beneke, P. Falgari, S. Klein and C. Schwinn, Nucl. Phys. B 855 
(2012) 695 [arXiv: 1109.1536]; S. Alekhin, A. Djouadi and S. Moch, Phys. Lett. B 716 (2012) 214 
[arXiv: 1207.0980]. 


29 


[11] ATLAS Collaboration, Eur. Phys. J. C 74 (2014) 10, 3109 [arXiv: 1406.5375]. 

[12] CMS Collaboration, Phys. Lett. B 728 (2014) 496 [Phys. Lett. B 728 (2014) 526] [arXiv: 1307.1907]. 

[13] A. J. Buras and M. K. Harlander, Adv. Ser. Direct. High Energy Phys. 10 (1992) 58. 

[14] M. K. Gaillard and B. W. Lee, Phys. Rev. D 10 (1974) 897. 

[15] J. R. Ellis, M. K. Gaillard, D. V. Nanopoulos and S. Rudaz, Nucl. Phys. B 131 (1977) 285 [Nucl. 
Phys. B 132 (1978) 541]. 

[16] Gutter Group Gollaboration, Eur. Phys. J. G 74 (2014) 3046 [arXiv: 1407.3792]. 

[17] M. Ciuchini, E. Eranco, S. Mishima and L. Silvestrini, JHEP 1308 (2013) 106 [arXiv: 1306.4644]; 
M. Ciuchini, E. Franco, S. Mishima, M. Pierini, L. Reina and L. Silvestrini, arXiv:1410.6940. 

[18] K. Seidel, F. Simon, M. Tesar and S. Poss, Eur. Phys. J. C 73 (2013) 8, 2530 [arXiv: 1303.3758]; 
T. Horiguchi, A. Ishikawa, T. Suehara, K. Fujii, Y. Sumino, Y. Kiyo and H. Yamamoto, 
[arXiv:1310.0563]. 

[19] M. Beneke, Y. Kiyo, P. Marquard, J. Piclum, A. Penin and M. Steinhauser, arXiv:1506.06864. 

[20] R. Barbieri, M. Beccaria, P. Ciafaloni, G. Gurci and A. Vicere, Phys. Lett. B 288 (1992) 95 [arXiv:hep- 
ph/9205238]. 

[21] J. M. Gornwall, D. N. Levin and G. Tiktopoulos, Phys. Rev. D 10 (1974) 1145; G. E. Vayonakis, 
Lett. Nuovo Gim. 17 (1976) 383; M. S. Chanowitz and M. K. Gaillard, Nucl. Phys. B 261 (1985) 379. 

[22] M. Misiak et ai, Phys. Rev. Lett. 114 (2015) 22, 221801 [arXiv:1503.01789]. 

[23] M. Czakon, P. Fiedler, T. Huber, M. Misiak, T. Schutzmeier and M. Steinhauser, JHEP 1504 (2015) 
168 [arXiv:1503.01791]. 

[24] B. A. Kniehl and M. Spira, Nucl. Phys. B 443 (1995) 37 [arXiv:hep-ph/9501392]. 

[25] D. Abbaneo et ai, EGFA High Luminosity LHC Experiments Workshop: Physics and Technology 
Challenges, ECFA/13/284. 

[26] L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 1945. 

[27] A. J. Buras, M. E. Lautenbacher and G. Ostermaier, Phys. Rev. D 50 (1994) 3433 [arXiv:hep- 
ph/9403384]. 

[28] Particle Data Group Gollaboration, Ghin. Phys. G 38 (2014) 090001. 

[29] H. Baer et ai, “The International Linear Collider Technical Design Report - Volume 2: Physics”, 
[arXiv:1306.6352]. 

[30] K. G. Ghetyrkin, J. H. Kuhn, A. Maier, P. Maierhofer, P. Marquard, M. Steinhauser and G. Sturm, 
Phys. Rev. D 80 (2009) 074010 [arXiv:0907.2110]. 

[31] Heavy Flavor Averaging Group (HFAG) Gollaboration, [arXiv:1412.7515]. 

[32] See G. Martinelli’s talk at the EPS 2015 conference and the UTfit website. 

[33] A. Andreazza et al, Frascati Phys. Ser. 60 (2015) 1. 

[34] T. Aushev et ai, arXiv:1002.5012. 

[35] LHGb Gollaboration, Eur. Phys. J. G 73 (2013) 2373 [arXiv:1208.3355]. 


30 


[36] A. J. Buras, D. Guadagnoli and G. Isidori, Phys. Lett. B 688 (2010) 309 [arXiv:1002.3612]. 

[37] A. J. Buras, M. Jamin and P. H. Weisz, Nucl. Phys. B 347 (1990) 491; 

[38] Gorbahn et al., work in progress. 

[39] J. Brod and M. Gorbahn, Phys. Rev. Lett. 108 (2012) 121801 [arXiv:1108.2036]. 

[40] J. Brod and M. Gorbahn, Phys. Rev. D 82 (2010) 094026 [arXiv: 1007.0684]; Nucl. Phys. B 523 (1998) 
40 [arXiv:hep-ph/9710245]. 

[41] J. Brod and M. Gorbahn, Phys. Rev. D 78 (2008) 034006 [arXiv:0805.4119]. 

[42] G. Isidori, F. Mescia and G. Smith, Nucl. Phys. B 718 (2005) 319 [arXiv:hep-ph/0503107]. 

[43] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68 (1996) 1125 [arXivdrep- 
ph/9512380]. 

[44] C. Lehner, E. Lunghi and A. Soni, [arXiviarXiv: 1508.01801]. 

[45] CMS and LHCb Collaborations, Nature 522 (2015) 68 [arXiviarXiv: 1411.4413]. 

[46] G. Buchalla and A. J. Buras, Nucl. Phys. B 548 (1999) 309 [arXivdiep-ph/9901288]. 

[47] C. Bobeth, M. Gorbahn, T. Hermann, M. Misiak, E. Stamou and M. Steinhauser, Phys. Rev. Lett. 
112 (2014) 101801 [arXiv:1311.0903]. 

[48] F. Mescia and C. Smith, Phys. Rev. D 76 (2007) 034017 [arXiv:0705.2025]. 

[49] J. Brod, M. Gorbahn and E. Stamou, Phys. Rev. D 83 (2011) 034030 [arXiv: 1009.0947]. 

[50] T. Blum, R. Zhou, T. Doi, M. Hayakawa, T. Izubuchi, S. Uno and N. Yamada, Phys. Rev. D 82 
(2010) 094508 [arXiv: 1006.1311]. 

[51] S. Borsanyi et al., Science 347 (2015) 1452 [arXiv: 1406.4088]. 

[52] N. Carrasco, V. Lubicz, G. Martinelli, C. T. Sachrajda, N. Tantalo, C. Tarantino and M. Testa, Phys. 
Rev. D 91 (2015) 7, 074506 [arXiv: 1502.00257]. 

[53] A. Alberti, P. Gambino, K. J. Healey and S. Nandi, Phys. Rev. Lett. 114 (2015) 6, 061802 
[arXiv:1411.6560]; T. Mannel, A. A. Pivovarov and D. Rosenthal, arXiv: 1506.08167. 

[54] RBC and UKQCD Collaborations, Phys. Rev. D 88 (2013) 014508 [arXiv: 1212.5931]. 

[55] G. Isidori, G. Martinelli and P. Turchetti, Phys. Lett. B 633 (2006) 75 [arXiv:hep-lat/0506026]. 

[56] G. Altarelli, R. Barbieri and S. Jadach, Nucl. Phys. B 369 (1992) 3 [Nucl. Phys. B 376 (1992) 444]; 
G. Altarelli, R. Barbieri and F. Caravaglios, Nucl. Phys. B 405 (1993) 3. 

[57] G. Degrassi, P. Gambino and P. P. Giardino, JHEP 1505 (2015) 154 [arXiv:1411.7040]. 


31 


